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DIRECT READING FLUID FLOUMETER 


by 


H.L. Cummings and &.E. Mallick 


Submitted for the degree of Master of Science in the 


Department of Aeronautical Engineering on May 18, 1951. 


Summary 


A flowmeter designed to directly measure mass flow 
as a function of the pressure difference across a rotating 
cylinder in a fluid stream was investigated. Calibration 
curves for the meter indicate an accuracy of measurement 
of two per cent. The useful range of the meter was found 
to be limited by the ratio of the cylinder tangential 


speed to the free stream velocity. 
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INVESTIGATION AND CALIBRATION CF A 


QIRECT BEADING FLUID FLOWYETER 


Ie INTRODUCTION 


A meter which vould directly indicate the mass 
rete of fluid flow in a duct would have many engineering 
arrnlications. The present standard flowneterine devices 
(the venturi meter, nozzle, and orifice) require the 
measurenent of three stream properties and the computa 
tion of the mass rate of flow by use of an equation and 
various calibration coefficients. The problem of the 
metering of a fluid would, therefore, be greatly simpli- 
fied if an indication of mass rate of fluid flow could be 
obtained from a direct resding instrument. Such a device 
would be of use in the control of fueleair ratios under 
conditions of varying density as experienced in altitude 
operations of aircraft and nissiles. In measuring the 
nass flow rate of eases of varying density or of unknowm 
counosition such 4 meter mirht be enrloyed. 

The rurrose of this thesis is, then, the investi- 


gation, calibration, and analysis of such a direct read- 











ing flovmeter. 


The flowneter upon which this renort is based was 
designed and constructed by Tennant and Turner and 
rerorted upon in their thesis 7 * at the Massachusetts 
Institute of Technolory in 1950. The original equirment 
has been modified by the present authors and this report 
includes results obtsined from the operation of the 
wal Gad! @rnanatue and analysis of these results. 

The flowmetor ciscussed in this report consists 
essentially of a square duct within which a cylinder ro- 
tates in a counter-clockwise dircction about an axis per- 
pendicular to the direction of flow as shown in Fig. l. 
There is, conseauently, a speeding up of the flow in the 
duct below the cylinder and a retardation of velocity 
above the cylinder. This difference in velocity results 
in a difference in static pressure at the walls of the 
duct directly above and below the rotating cylinder. 
The mapnitude of this pressure difference is directly 
related to the mass flow in the cuct as will be further 
discussed in Chapter II. 

The principle upon which this device is based, that 


of a rotating object in a fluid stream, was originally 





* Numbered superscripts refer to the Bibliosraphy, pe 
LO and AL . 














discussed by Marnus 2 in 1851 in connection with the 
flight of projectiles and for this reason is generally 
known as the "hapnus sffect™. The problem was also 
further discussed by Lord Rayleifh 3 in 1877 and by 
Lafay ? in Paris in 1912. In the early 1920's the 
rotating cylinder was considered by many aerodynamicists 
as a more pronising liftine device than the presently 
accented and proven Joukowski airfoil. This outlook was 
based primarily on the potential flow analysis of the 
problem which theoretically predicts a lift coefficient 
for a rotating cylinder of #77 which is considerably 
higher than that obtained fron airfoils of the Joukowski 
type. An explanation of this analysis is piven by Gold- 
stein * , 

hxtensive experimentation was carried out at the 
German laboratories st Gottingen under the direction of 


Prandtl ime ‘ 


Flow pictures obtained by Prandtl are 
Shown in Fig. 17 and are ciscussed in cetail in later 
names of this rerort. 

The promisine results obtained at Gottingen led 
Flettner to anply the vrinciple of the "™vapnus Effect" 
to ship propulsion; replacing the sails with two 45 ft. 
high, 9 ft. in dianeter, rotating cylinders. The 


Flettner rotor ship made several semi-successful mins; 


== 
i; dl oe 
ss > 
= =p 








however, the auestion of the practicability of such a 


mode of propulsion was a source of much controversy. Two 


papers by Betz 11 and Ahlborn wi discuss the pros ane 


cons of the Flettner rotor shin. 
Concurrent with the Gottingen experimentation, 
considerable research concerning the rotating cylinder 


in a real fluid was carried out by Thom 6,7,8 


in Fneland. 
vith the failure of the Flettner rotor ship to 


fulfill the exnectations of its desifners, a con- 





secuent lacx of interest in the rotating cylinder eas a 
lifting device resultéd. In 1925 Reid 9 of the NACA 
described experizents with a rotating cvlinder in an 
airstream which confirmed much of the data obtained by 
other investigators. 

To date the practical applications of the lifting 
force obtained fron the "Mapnus Fffect" have been larrely 
unsuccessful. Ferhaps the information presented in this 
report will result in a practical enployment of this 


long knovm and nuch studied nhenomenon. 











II. THEORETIC ANALYSIS 


2e1 Introduction. The develoment of the theoretical 
flowmeter equation contained in this chapter follows 
substantially the original derivation as set forth by 
Tennant and Turner a This treatment is included in 
this report, however, both for completeness and for 
correction of a few minor errors, notably an error in 
the calculation of the nunerical constant in the 
theoretical flowmeter equation of the original thesis. 
The method followed in the derivation is essentially an 
extension and combination of two analyses as reported 
by Durand 10, 

The meter is basically as shorm in Fir. 1. By 
neesuring the pressure difference AF between points 
A and B when the cylinder is rotating at a knowm 
constant sneed it is proposed to heasure thé mas rate 
of flow through the apparatus. For realization of such 


a meter there must exist a relation of the form 


WV = Const. x AF (1) 


where Via denotes the mass rate of flow. Expressing Mz 


as the product of f , density; A, area; and V, velocity 


we obtain PAV= Cohst. « AP 











which, with a known constant meter peonetry, may be 


expressed as 


AP=const x PV. (2) 


Thus it remains to cenonstrate the existence of relation 
(2) and to obtain a nuverical value for the conmuemt of 
relation (1). 

For the accomplishrient of the above ains several 

Simplifications of the actual flow are assuned: 

a. The flow is two dimensional. 

b. The flow is irrotational. 

c. The fluid is a non-viscous, incompressible, 
ideal fluid. 

d. The rotating cylinder imparts to the fluid a 
circulation | , equal to 27T tines the 
product of the radius, R, of the cylinder, 
anc the peripheral speed, U, of the cylinder. 

The above assumptions, combined with the require- 

ment that the flow be continuous, make possible the 
employnent of the theory of potential flow for synthesis 
of the flow pattern within the eter. Having an express- 
ion which describes the flow in the eter, it is then 


possible to find the velocities at points A and B. With 











these velocities known, the Bernoulli equation is applied 
between the tio pnoints resulting in an expression relating 
the difference in pressure between A and B to the stream 


velocity and density. ; 


2.2 Kethod of Analysis. The potential flow pattern 

about a stationary cylinder in a fluid stream is produced 
by combination of a sourceesink doublet and parallel flow 
(Fig. 2a). If a potential vortex, the core of which is the 
cylinder (Fie. 2b), is added to this picture, the flow 
pattern for a rotating cylinder in an unrestrained fluid 
stream is the result (Fig. 2c). In the actual meter the 
flow is restrained by the walls of the meter duct. The 
effect of the walls upon the flow pattern is simulated by 
the method of images originally devised by Glauert 26 | The 
wall correction is avplied separately to the flow patterns 
of the stctionary cylinder (Fir. 3a) and the potential 
-vortex (Fig. 3b). The supervosition of these two correcte 
ec patterns is then the true potential flow pattern for 


the meter. 


23 The Stationary Cylinder and Its Images. Consicer 


the courlex 2 plane with x and iy axes as showm in Fig. 
3a. The conrlex variable 2-=x + iy describes any point 


in this nlane. The stationary cylincer is located with 





8. 


center at the origin. The effects of the duct walls on 
the flow are simulated by an infinite series of cylinder 
images extending in both positive and negative directe 
ions. Centers of the image cylinders are located on the 
jy axis at iy = ¥ ink where n is a positive integer 
and k is the duct height. 

The potential funetion / = Drel completely 
describes any potential flow. p ' the velocity poten- 
tial, Y ; the stream function, and F are functions of 
the complex variable 2. The potential function for the 
flow about a stationary cylinder in a parallel stream 
of velocity V is * 


fk =Vz + (3) 


where C is a constant. This function describes the flow 
at any point z in the complex plane due to the obstruct- 
ing cylinder at the origin. Similarly 


- C 
fF = lz * Funk 





* Hildebrand, Ref. 13, pe 547, Eqe 17h- 











describes the flow at any noint 2 due to 2 cylinder 


image at iy = + ink and 


_ C 
fF = Vz + ZHiNMNK 


describes the flow at z sue to a cylinder innge at 


ly = — ink. The potential function for the collectivity 


of images and real cylinder is then 


= Vz +C[#+ > (eime + sine | 


K 
77 


F=VerE {catey + 0E)> [ariel 


Making the substitution C = c 





and sinmnlifving: 


The sun of which, by Fittag-Leffler's theorem *, is 


fF = Vz +C, coth = (5) 





* Jeffreys and Jeffreys, Ref. li,, pe 357, Ee 8 





The stream velocity wat z is given by the 


dF 
dz 


derivative at that point 


pg 
_ ot ee 
Sink? Te 


To evaluate a » consicer the physical condition 
that the velocity is zero at the points on the x axis 
where 2 =f R, R being the radius of the cylinder. 

Then 





é, = £ sinh® BA 


The velocity at any rnoint 2 is then 


R 
W = V// _ sie = 


2 72 
SIMA —— 


a | 
At the walls ‘ireectly above and below the cylinder where 


jie wee , the velocity is 











27R 
SIAh We 


w=V /- 
sink ti 


which reduces to 





2 Fe (8) 


WAV + SIN 


which is the velocity at points A and B on the walls, 
directly above and below the cylinder, due to the 
obstruction to the flow by the cylinder. To this must be 


aaced the velocity due to the rotation of the cylinder. 


204 The Vortex _and Vortex Images. Again consider the 


comnlex vlene as shovm in Fig. 3b. The circulation 
suprlied to the fluid by the rotation of the cylinder is 
represented by a rotential vortex at the origin, the 

core of which is the cylinder. The effect of the duct 
walls is sinulated by ean infinite series of vortex images 
located on the iy axis with centers at ly =~ ink. To 
give the correct velocity distribution transverse to the 


flow direction, alternate vortices are piven opposite 








rotations; therefore, there is a sistance 2k between 


vortices of the same circulation. 


The votential Punction at any roint 2 for a vorte@r 


located at the oririn is * 


—_ 


where / °, the circulation, is vositive in the counter- 


clockwise direction. 


It will be mathematically convenient if we sun 





velocity expressions, therefore, taxing the derivative 


dF 


= _—o- = TL (10) 
VA/ dz 


— 
CML. = iE 


which is the velocity at any point z cue to a vortex of 
Strengt* | locatec at the orifin. 
The velocity at any noint 2 due to a vortex of 


Bane diréction of rotation locatéd at iy © # inzk is 


.c) ae 
a) 5 bal ee 0/7) 





® Wild@ebrand, Ref. 13, pe 565, Uae 172. 








The velocity at any point @ due to a vortex of 


same rotation at iy = ein2k is 


ja y 
CITL Z #F4h2K 


Ay = 


The sum of the velocities at % due to the real 
vortex and a vositive and negative infinite series of 


vortices of the same direction of rotetion is 


me / y / 
AE 2 jad ™ 2 Zhen . = Lien 





Making the substitution Ve = / AK » and simplifying: 


Eg +0) [acre el} 


The sum of which by IiittageLeffler's theorem * is 


/ Oe (12) 
: — / 
W = coth ae 





* Jeffreys and Jeffreys, Ref. 14, pe 357, Fa. 8. 











Revriting in terms of / 3: 


F#KRL 2K 


which is the velocity at any point 2 due to the real 
vortex and the infinite series of images of same directe 
ion of rotation. 

To this itust be added the velocity due to the 
images of oprosite rotation. 

The velocity at 2 due to a vortex of opposite 


rotation at iy = #L(2n-IK is 


y} 
Z-L(2n-1) K 





The velocity at 2 due to a vortex of opposite 


rotation at iy = 22s K is 


=f / 
CITA ZF LeA-YK 


id 


AZ 


The sun of the velocities at z due to the posite 
ive and negative infinite series of vortices of opposite 


rotation is 








~~ 





Be) 
= FET > ae Z-< (2n-)K ep. Se / 


Making the substitution y a; 4K, anc sinplifying: 


4 7 Z T VA 
—_ ra 2) > |e carnal 
/ 2k | 


The sun of which, again by littag-Leffler's theorem *, is 


wW Ta, Ve 
W L ca LE 
Rewriting in terns of / 


= 6 
W/= a tanh Le (16) 


which is the velocity at 2 cue to the infinite series of 
vortices of onrosite rotation. 


Combining, eountions (13) and (16) the total ve- 


* Jeffrevs and Jeffreys, Ref. 14, pe 357, ta. ll. 











locity at any point z due to the real vortex and all its 


images is 


ae WZ WZ 


which can be simplified to 


a i / 
aeeey Sthh WE 
K 





At point A on the wall directly above the cylinder 
_K 
Ze tA and the velocity due to the vortex and wall 
effect upon the vortex is 


fr 


= — 24 (19) 
_ 


Similarly at point B the velocity is 


lam 
= += 20 
LA/ 5 (20) 


To these velocities must be added the velocity due 
to the stationary cylinder and its images. Combining 
equations (8) and (19) and replacing [| by27RU we 











obtain for the total velocity at A: 





vas 
W, = V//4 Sinha cs U (23) 


Combining equations (8) and (20) and similarly 


replacing we obtain for the total velocity et Bs: 


We =V(/ + Sinh + FRY hie 


2e5 Theoretical Flowmeter Equation. With the velocities 
at points A and B known, we apply the Bernoulli equation 


between points A and B. 


Py 2 
aP=Z-(wi-w = 


which results in 


SF 477R 
= (24) 
AP=s9 WV (+ Sihh 
Equation (24) vith a constant speed of cylinder rotation 


is of the forn 


AP =const. xPV (2) 








which establishes the theoretical basis for the flowmeter. 
Combining equation (24) with the usual equation 
for mass flow and substituting the dimensions of the 


meter we arrive at the theoretical flowmeter equation 


dees NT (25) 
Ma = ./55 x lO Yyy— 


in which M, is in lbm/sec, Af is in inches of water, 


anc N is the revolutions per minute of the cvlinder. 














Jif. &  ULPMENT AND PROCEDURE 


The flowneter investigated consisted of an alumi- 
num cylinder four inches in ciameter mounted in an eirht 
inch square duct as shown schematically in Fig. l. 
Annular aluninunm end plates were inserted in the sides 
of the duct walls to house the ends of the cylinder. This 
arrangement of end plates was provided to eliminate the 
flow of air around the ends of the cylinder; thus sirme- 
lating two dimensional flow. The shaft ends of the cyline 
der were fitted into New Departure, S.K.F. No. 2, self- 
aligning ball bearings. In order to assure accurate 
alignment of the cylinder within the duct walls, the 
races of the ball bearings were housed in four inch 
square, 1/4 inch thick, aluminun plates that were attach- 
ed to the duct walls. The photorraph in Fir. 4 shows the 
details of the cylinder mounting in a cut away view of 
the test section. Since sealed ball bearings were not 
used, sore leakage of air was expected around the cyvlin~ 
der shaft. During the test runs, however, the pressure 
in the test section wis held as closely es possible to 
atmospheric in order to mininize leakage. 

The cylinder shaft was extended through the test 


section w2ll on the power input side, and courled to the 














rower dirive by -eans of a flexible coupling. The photoe 


graph in Fig. 5 shows the details of the power drive 
installation. 

A standard Dunmore No. 5, 1/2 hep. grinding motor 
was used to rotate the cylinder. The quill of the gina 
ing .otor was coupled to the cylinder shaft as show in 
Fig. 5. The Punore No. 5 grinding motor was chosen to 
drive the cylinder because it afforded a suitable range 
of rotational speeds (up to 14,000 RPM for the equipment 
used) for testing purposes. Since accurate alignment of 
the driving motor quill with the cylinder shaft was 
necessary in order to eliminate vibration problems, the 
“sotor was mounted on a bracket that afforded easy adjuste 
ment for alignment purroses. The motor mounting bracket 
permitted adjustment in two directions at right angles to 
each other in the horizontal plane. A vertical adjustment 
on the driving motor permitted alignment in the vertical 
direction. Positive control of the RPM of the motor was 
obtained by controlling the voltage input to the motor 
with a Variac. The RPM of the cylinder was measured with 
a Strobotace 

Air for the flowneter was supplied in a closed 
circuit system by the DeLaval Air Compressor of the 


MeI.eT. Gas Turbine Laboratory. This compressor is normally 














used for the Gas Turbine Laboratory variable density 


wind tunnel anc is rated at 15,500 ft. /mins. Supply and 
return of air to the flowreter was afforded bv neans of 
a@ ten and one-fourth inch spirally welded pipo. A plan 
view of _———e is show in Fiz. 6, and a photogranh 
of the installation is showm in [Fiec. 7. 

In conducting the tests it was found to be nore 
practicable to run the cylinder at a constant RPM, and 
to vary the air flow. Seven different RPMs were selected 
for making the test runs: 3,300; 4,350; 6,400; 8,700; 
10,000; 12,000; and 14,900. The selection of these RPMs 
was arbitrary; they were chosen to illustrate the effect 
of RIM on flow measurenent for this type of flowmeter. 
The derendence of the flowmeter on the ratio of cylinder 
tancential speed to free stream velocity for accurate 
flow measurenent is discussed in Chapter IV. 

Calibration of the flowneter was effected by use 
of a standard ASME square edged orifice which was placed 
in the flow circuit as showm in Fig. 6. The diameter 
ratio of this orifice was 0.6. Calibration curves for the 
orifice are shown in Anpendix II. The pressure difference 
across the orifice was measured with a water manometer at 
low air flows with a means provided for switching to a 


mercury manometer at high air flows. 











Screens and brass flowestraightening tubes of two 


inch outside diameter were placed in the flow circuit at 
two places (see Fir. 6) so that parallel flow conditions 
would be more closely arproximated upstream of the flowe 
meter anc the orifice. 

The mass flow of air was varied fron approximately 
0.6 1b./sec. to 6 1b./sec. for each test run. For the 
equipaent used, this mass flow of air corresponds to free 
strean velocities varying from aprroximately 20 ft./sec. 
to 190 ft./sec.. The hirhest free stream velocity corres- 
ponds to a Mach Number of approximately .17. It is assumed, 
therefore, that conpressibility effects were negligible. 
The flow Reynolds Nunber (based on the dianeter of the 
orifice) varied from .14 x 10° Co 1.2 xX 10° for the range 
of free stream velocities used for these tests. 

Two different locations of the static pressure 
tans for measuring the pressure difference across tre 
rotating cylinder of the flowmeter were investig:rted. 

One location was in the duct walls directly above and 


below the center of the cylinder as shown schematically 





in Fig. 1 as points A and B. The authors were guided in 


their assumption that this was the optimun location * for 





* Betz, Ref. rls. Pe 18. 











the taps since Prandtl reported in his experiments with 


rotating cylinders that theoretical potential flow 
conditions were more closely approximated in this region. 
The pressure taps located in this position are referred 
to in this report as the primary pressure taps. The 
second location investigated was at the end of the cylin- 
der with the taps diametrically opposed at the surface 
of the cylinder. The pressure taps located in this posite 
ion are referred to as the secondary pressure taps. The 
location of these taps is showm schematically in Pig. 1 
at points C and D. 

The pressure difference from the primary taps 
was read on an Ellison Inclined Vertical Draft Gauge 
which was calibrated to read pressure difference directly 
in inches of water. Since the range of this instrument 
extended only to 12 inches of water provisions were made 
for switching to a conventional water manometer when 
larger readings were necessary. A photograph of the 
manometer board used for taking readings is shown in 
Fig. 3. 

Fressure difference from the secondary taps was 
read on a conventional water manometer. 

In order to provice a neans for evaluating static 


pressure lozs across the flowmeter, static pressure taps 





2h, 


were installed 12 inches upstream and 12 inches down- 
stream from the rotating cylinder as shown in Fig. 6. 
The difference in pressure was measured on a convention- 
al water manometer. 

in the test procedure the cylinder was rotated at 
a selected RPM. This RPM was kept as constant as possible 
throughout the run by using a Variac to control as 
necessary the voltage input to the motor. The mass flow 
of air through the system was started at as low a value 
as could be obtained with the test equipment and ine 
creased in selected increments. Simultaneous readings 
of all manometers were taken during each increment of 
air flow with sufficient tine being allowed between 
readings for steady flow conditions to be established. 
The mass flow for each flow increment was computed from 
the orifice data, and the reading of the flowmeter was 
calibrated with this mass flow. The data taken and 


sample calculations are shown in Appendix I. 











IV. DISCUSSION OF RESULTS 


4.1 Introduction. A graphical dnalysis of the results 
of the investigation and calibration of the flownieter is 
given in Figs. 9 through 16. Firs. 9, 10, and 1l present 
curves relatinz actual mass fluic flow anc meter reade 
ings. These curves incicate the true worth of the apna- 
ratus. Figs. 12, 13, and 14 present plots involving 
various parameters and non-dinensional quantities. These 
curves are offered in partial explanation of the results 
portrayed in the first groun of plots. These curves also 
show the linitations of the meter anc offer guidance in 
selecting the meter characteristics for an expected range 
of mass flow.Fies. 15 an’ 16 present other data obtained 


in the course of this investircation. 


l.2 The Velocity Ratio U/V. The non-dimensional velocity 
ratio U/V is a quantity which figures extensively in 
correlation and explanation of mcter cat. In this ratio, 
U denotes the peripheral velocity of the rotating cylin- 
der, while V cenot:s the free undisturbed stream velocity 
ahead of the meter. That the factor U/V exerts a rreat 
influence on the flow picture within the meter is denone- 
strated both fro: potential flow theory and fron consider- 


ation of real fluid flow phenomena. The potential flow 











pattern for a stationary cylinder in parallel flow as 


shown in Fig. 2a is, in effect, the flow pattern for a 
U/V of zero. The flow stagnation points are seen to 

exist exactly at the front and rear of the cylinder. If 

a circulation is added, U/V becomes a positive quantity 
and the stagnation points move along the surface of the 
cylinder closer to the region of decreased velocity as is 
pictured in Fig. 2c. Further increases in circulation, 
that is, further increases in U/V, result in junction of 
the two stagnation points on the underside of the cyline 
der and movement of the stagnation point into the fluid. 
The effect of the ratio U/V on the flow pattern of a real 
fluid can be clearly seen in Fig. 17 which is a series of 
flow pictures taken by Prandtl at Gottingen. For values 
of U/V¥ at which the flow picture of the real viscous 
fluid resembles the potential flow pattern, correct in- 
dications by the flowmeter should be expected. when, how- 
ever, the viscous effects are such as to cause conplete 
disagreement between the actual pattern anc the pattern 
of the potential flow, we should expect that our meter 
reading will not be an indication of the actual mass flow. 
Examination of Fig. 17 would lead us to expect no worth- 
while meter indications below a U/V of about one. 


This critical value of U/V cannot be applied di- 











rectly to our flowmeter, however, cue to the difference 


in measurement of the quantity V. The frandtl pictures 
are of a cylinder rotating in an unrestrained fluid 
stream. The flowmeter stream is considerably restrained 
by the duct walls; consecuently, the true VY is greater 
than the free stream V used in conruting the U/V of this 
rerorte Thom? gives a method for correction of U/V to 
include duct wall influence; however, the authors doubt 
the validity of the approximations made and have omitted 
any corrections to U/V in this paper. Thus, the actual 
values of U/V computed in this rerort may not be compared 
with values obtained by other investisators, but U/V 
still remains an imnortant parameter for correlation of 
flowmeter data at various cylinder rotational speeds and 
mass flow ratcs. 

That the ratio U/V is a determining factor in flow- 
meter performance may be deduced in another manner. In 
problems of flow sinilitude, the Reynolds Number, which 
is the ratio of inertia forces to viscous forces, has 
been theoretically and experimentally shovm to be aleone 
trolling parameter. The flow in the meter may be consider-~ 
ed as the superposition of two flows, the parallel flow 
in the duct and the circulatory flow. The Reynolds Number 


for the duct flow may be written as 








PVe 





a 


where /7 is the fluid density, YV is the stream velocity, 
_“©& is the fluid viscosity, and ¢ is the hydraulic 
diameter of the ducte The Reynolds — for the circu- 
latory flow may be written as 
R,=-LU¢ 

5 eal eee 
where U is the peripheral velocity of the cylinder and d 
is the cylinder diameter. Thus, for the sane fluid and 
fixed meter geometry, the factor U/V is seen to be a 
ratio of Reynolds Numbers and hence viscous effects 


should be expected to be a function of U/V¥. 


4.3 Discussion of Figs. 9, 10, and il. Figs. 9, 10, and 


ll are plots of actual mass flow, Mae, versus a which 


is the meter reading in inches of water divided by the 
A 


Z 
cylinder RP. The abscissa is actually 7,7 x /O to 
enable the use of whole nunber scale divisions in plote- 
ting. The theoretical flowneter equation derived in 


Chapter II] 


ee 
Mig = SS SK NO Fy (25) 
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is plotted on all three graphs as a dashed straight line. 

Fig. 9 shows the performance of the meter as 
compared with the theoretical prediction for 5 RPMs, the 
highest of which is 10,000 RPM. The useful range for any 
RPM is represented by the nearly straight line portion of 
the curve extending outward from the origin. It can be 
observed that, at low RPMs, this useful range is extreme-~ 
ly narrow and widens with increasing RPM. The curve for 
10,000 RPM shows a range of useful operation up to a masa 
flow of about 3.9 lbm/sec . 

Fig. 10 is a similar plot which shows the useful 
range of operation for an RPM of 12,000. Three independ- 
ent runs made on different days are plotted on this graph 
to demonstrate the reproducibility of data. The three 
runs may be represented by a mean line up to a mass flow 
of about 4.6. This mass flow corresponds to a U/V of 
about 1.6 which will be shown in the discussion of the 
next section to be the critical value of U/V¥ below which 
the meter is not accurate. Therefore, divergence of the 
data at this point as occurs on this graph is to be exe 
pected. Examination of the data for the three runs, for 
mass flows above the critical value of U/V, shows a 
maximum deviation of 2% from the mean line. 


In Fig. 11 three separate runs at 14,000 RPM are 
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plotted. A maximum deviation of about 2% from the mean 
line is also found at this RFM. 

Both Figs. 10 and ll show decided depressions in 
the curves between mass flows of 1 and 2 lbm/sec . These 
depressions occur at high values of U/V of the order of 
6 to 9. It is suspected that in this region a component 
of the flow velocity is directed into one or both of the 
static pressure taps, thus causing a slight irregularity 
in the calibration curves. 

Figse 9, 10, and ll are the calibration curves for 
the flowmeter which might be used in actual operation; 
however, it is necessary to examine other plotted re- 
lations for further insight into the meter flow process, 
4ey Discussion of Fig. 12. The velocity ratio U/V is 
oa x/O7 4n Fig. 12 for the seven RPMs 





plotted versus 
which have been studied. Also plotted is the theoretically 
predicted line for 4,350 RPM. The theoretical curves for 
the other six RPMs have been omitted since they would 
conflict with the curves plotted from the actual data. 

The theoretical curves are a family of hyperbolas all 
similar to the curve for 4,350 RPM. According to the 
theoretical relation derived in Chapter II, sf is di- 


rectly proportional to Wa , the mass rate of flow. With 
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a constant fluid density and duct area the velocity V is 
also directly proportional to the mass flow. For a constant 


RPM, U, the peripheral velocity, is a constant and U/V is 





M AP 
thus proportional to Baa wy is proportional to “a2 ; 


therefore, we are, in effect, plotting a relation of the 
type "x versus 1/x" which is, of course, a hyperbola. 

The actual curves are displaced from the theoreti- 
cal curves due to the many simplifications of the flow 
made in the theoretical calculation. The curves for 10,000, 
12,000, and 14,000 RPM very closely conform to the hyper~ 
bolic shape which is again an indication of their useful 
range of operation. The significant feature of Fig. 12 
is that the curves for all RPM below 12,000 show a break 
at a U/V of about 1.6. This confirms the suspicion that 
the flow process within the meter is primarily dependent 
upon the value of the velocity ratio. Above the critical 
value of U/V the speed of rotation of the cylinder is 
sufficient to prevent wapiranhes of the boundary layer 
fluid from the cylinder and consequent formation of 
vortices. Here, then, the flow pattern still resembles 
the potential flow pattern. At the critical value of U/V 
the rotation of the cylinder is no longer sufficient to 
prevent boundary layer separation in the rear of the 


cylinder anc consequently the flow picture is no longer 
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similar to that predicted by potential flow. Probably in 
this region of operation the static pressure taps at "A" 
and "B" no longer measure only static pressure but now 
recover a component of the velocity directed toward them 
in the turbulent and presently unpredictable flow which 
results when separation from the cylinder occurs. The 
flow pattern at and below this critical point is probably 
much like the photographs for U/V 's of 1/2 or 1 of Fig. 
17. Stagnation of a component of the velocity in the 
pressure taps below the critical value of U/V is probably 
also the cause of the change in sign of the measured 
pressure difference which occurs for 3,300 RAM and which 


is shown on both Figs. 9 and 12. 


405 Discussion of Fig. 13. Fig. 13 is a plot of cali- 
bration coefficient, K, versus U/V and is included here 
for academic consideration rather than for actual use in 
connection with flowmeter operation. For actual flowe 
meter operation, use of a calibration curve such as Fig. 
11 is recommended. The calibration coefficient is de- 
fined as the ratio of the actual mass flow to the theoe 
retical mass flow for a given value of = « Thus, it 


is a ratio of the instantaneous slope of the actual 


curves of Figs. 9, 10, and ll to the slope of the theoe 
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retical flowmeter equation. 
Below a U/V of 1.6 the values of K are meaning~ 


less and in many cases off the scale of this graph due to 
the predominance of viscous effects which bring about 
flow separation as discussed in the previous section. 
Above a U/V of 1.6 the data for various RPMs can not be 
said to plot as a single line; however, the general 
Shape of the curves for all the RPMs is similar. At high 
U/V 's there appears to be more correlation of the data 
for different values of RPM. Above the critical value of 
U/V the calibration coefficient is of the order of 4. 
This means that the pressure difference obtained for a 
given mass flow and RPM is only one fourth of the theoe- 
retically predicted pressure difference. | 
This result is in close agreement with data ob- 
tained by Thom 7 who measured the circulation produced 
by a rotating cylinder at a U/V of 2, The results of this 
experiment are shown in Fig. 18 in which circulation is 
plotted versus radial distance from the surface of the 
cylinder. Thom found that the actual circulation imparted 
to the fluid at a small distance from the cylinder is 
equal to about .24 of the theoretical value and remains 
constant at this value further from the cylinder. Since 


the circulation (/°=277RU ) appears as a factor in 








the theoretical equation for OF (Eq. 24, Chap. II), 


this would mean an expected reduction in MP to about 
one fourth of the original theoretical value. A value of 
K of about four is then to be expected from consideration 
of Thom's results. 

Exanination of Fig. 13 shows an increasing value 
of K at the hicher values of U/V¥. This indicates that 
the actual pressure difference measured by the meter is 
less than the predicted pressure difference by an even 
greater anount then the condition discussed in tte 
preceding paragraph. This further lowering of OP could 
possibly be due to a tilting of the flow pattern in the 
direction of flow. \hen the flow pattern is tilted, the 
meter pressure taps no longer measure the maximum prege 
sure difference and consequently a higher value of calie 
bration coefficient will be required. Tilting of the flow 
pattern at high values of U/V might also result in a 
component of the velocity into the pressure taps as 


previously discussed. 


4.6 Discussion of Fig. 14. The non-dimensional quan- 

tity = is plotted versus U/V in Fig. 14. q is the 
2. 

free stream dynamic pressure ( + V ). Below a U/V 


of 1.5 there is no correlation of data for various RPMs as 








expected; however, above the critical U/V, all the data 


for various RPMs appears in general to conform to one 
line. The factor ry is similar to tke lift coefficient 
used to describe airfoil characteristics. Thom 6 plots 
lift coefficient versus U/V for a rotating cylinder and 
obtains a curve very similar in shape to that of Fig. lk. 
Similar curves are also reported by Goldstein ‘“ from data 
of other investigators. The factor U/V is actually a 
ratio of Reynolds Numbers as previously discussed. This 
curve is then not constructed for any specific flow 
Reynolds Number but includes many varied values of this 
parameter. From this curve it may be concluded that flow 
in the meter is not a function of the flow Reynolds 
Number but is dependent only upon the value of the 


velocity ratio U/V. 


4.7 Discussion of Fig, 15. To investigate the optimum 


location of the meter pressure tars, secondary tars were 
installed close to the surface of the cylinder at one 
end of the cylinder ( points "C"™ and "D" of Fig. 1). 
Curves of Me versus as for these taps were meaning-~ 
less and are not included in this report. Fig. 15 

shows U/V versus <— for these secondary pressure tapSe 


The shape of the curves for the various RPMs is seen to 








vary with U/V as before. Breaks in the curves appear at 


constant values of U/V which is another indication of 


the dependence of the flow upon this parameter. 


4.8 Discussion of Fig. 16. The difference in static 


pressure ahead of and behind the flowmeter (12 inches in 
each direction) is plotted versus mass rate of flow in 
Fig. 16 for an RPM of 14,000. A curve of static pressure 
loss versus mass flow rate for a square edge orifice 

with pipe taps of area equivalent to the flowmeter 

passage area is presented for comparison. Flowmeter static 
pressure loss is seen to increase rapidly with mass flow; 
however, the magnitude of this pressure loss, even at 
high mass flows, is not exceptionally great, and is 
considerably less than that for the equivalent area 


orifice. 











Ve. CONCLUSIONS ARD RECOMMENDATIONS 


5el Conclusions. The practicability of a direct read- 
ing fluid flowmeter has been demonstrated anc calibration 
curves for various values of RPli have been presented. A 
maximus: error of 2% has been shown to exist in the results 
of three inderendent runs mace for each of two values of 
RPE. The derendence of the flowmeter upon the velocity 
ratio U/V has been discussed and various plots presented 
for substantiation of this discussion. It is concluded, 
therefore, that a flowneter of this type is capable of 
accurate and poncnelt measurenent of fluid flow within 
the discussed limitations. Due to the critical nature of 
the flow within the meter it would be necessary, however, 
to calibrate each individual meter before practical 


employment. 


5-2 Recommendations For Future Study. The tests mace 


in this investisation were restricted by the linitations 
of the eauiment used. In the theoretical derivation it 

has been shorn that the mass flow reading is inderendent 
of the density of the fluid. The wooden construction of 

the flovmeter prevented pressurization and consenuent 


density variation; thus, this theoretical property of 
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the meter has not been proven experimentally. 

The question of power consumed by the rotation 
of the cylinder is another avenue not completely investie- 
gated. The meter required the total output of the 1/2 
horsepower motor at 14,000 RPM for the maximum rate of 
mass flow studied. Detailed measurements of the torque, 
power required, etc. at various RPM and mass flows were 
not made. 

The problem of head loss attributable to the meter 
has only partially been investigated. Velocity traverses 
ahead and after the flowmeter combined with the static 
pressure readings at these points would probably present 
a more complete aspect of the losses (if any) connected 
with a meter of this type. 

The subject of optimum flowmeter dimensions is 
another possibility worthy of future investigation. In 
the meter studied the ratio of duct height to cylinder 
diameter was two. A greater value of this ratio would 
result in less obstruction to flow and less duct wall 
interference; however, restricting the duct size to 
practical limitations, this would mean a smaller cylin- 
der diameter and require higher values of RPM to maintain 


the proper range of the velocity ratio U/V. 





There is, then, a broad field of investigation 
ahead for perfection and complete explanation of this 


type of equipment. 
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GhOSSAP OF* SYi°BQLS 


pressure difference 

air mass flow (lb./sec.) 

air density 

area 

air stream velocity (ft./sec.) 
circulation 

cylinder radius 

eylinder tanrential sneed (ft./sec.) 

horizontal axis, imaginary vlane 

vertical axis, imaginary nrlane 
positive interfer 

duct heicht 

potential function 
comnlex variable 
constant 

strean velocity = dF/dz 
gravitational constant 
RPii (revolutions per minute) 
Reynolds Number 
cylinder ciameter 

fluid viscosity 


hydraulic diameter 








X 


cynanic pressure (1/2pv*) 1b. /ft2 
( Ma AcTVvAal 


calibration coefficient 
a THEO 
tenpsrature (cegrees Rankine) 
pressure 
orifice diameter 
pipe dianeter 
pounds mass 
sec seconds 
prinary refers to pressure tap location (see Fig. 1) 


secondary refers to pressure tap location (see Fig. 1) 














APPENDIX I 


SA SPL CALCULATIONS 


The sanvle calculations which follow are for Line 1, 





Table 1, vases 47 and 48. 

Column (1): Py orifice 
Py orifice = Py gauge + F 
P, gauge from manometer = O (not shown in table) 


Potm = 29085 in. hg. = 407 in. water 


Column (2): OP orifice 


OP orifice from manometer = .9l1 in. water 


Colwan (3): T 
Stream temperature used in computation of 


mass flow = 520 R. 


Govan (4): My oesmal 
Mass flow computed fron orifice data 
ae From Fig. Il-B, Appendix II, read K= .656 
for4OP orifice = .91 in. water 
be From Fig. II-A, Appendix II, read Y = 1.00 


for 4Pp orifice — .91 — .0022 
Pr orifice jt 








ce Compute MN. from: 


e 
Who >. H4#S DKY. = BO Poetrice 
where D5 = orifice diameter - 6.30 in. 
Py 29.85 in. hg. 
fT = 520 °R 


AP ar ee ~ 91 in. water 


m, = ( 1145(6.3)°(656) (10) pod (VW) = Ge 


Column (5): V 
zs ~ Ma _ 68 = 


= flowmeter duct area = ——— * a 


= 
’ 


Column (6): 4a 


q + / a > (002378) (19 25) “= FTF 


Column (7)! A Ppyowmecter (primary) 
from draft gauge MP = .32 in. water 


Column (8) N 
N = RPM measured by strobotae = 3290 RPM 





* W.A. Leary, Ref. 15, pe 9- 

















Column 


Column 


Column 


Column 


Column 


Column 


Colurnn 


(9): 


(10): 


(11): 


(13): 


(14): 


(15): 


U 
_-2HRN _ am (YeN329o) _ STF 


60 60 
where R= radius of cylinder = 1/6 ft. 


AP x 10* (primary) 


4 
¥ 


32 x10*= .974 


U/V 

U/V = olin = 2,88 

Op 

gq 
Ap = .32 in. water = 1.57 lbs/ft* 
AP 5967 = 3052 

q of fh, 

K 

K = calibration coefficient = We Aaa 

Ma THEO 
M, actual ='.68 lbm/sec 
b 

V, theo. = (.155) (qe x 10") = (.155)(.97h) 
K = 68 = eau 

AP flowmeter (secondary) 


from manoneterJOP = .51 in. water 


AP x 10” (secondary) 
‘al 


ae 1Ot= Ieee 
3200 
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(1) 


F? ABS 


ORIFICE 


“fH, O 


1,07 
1,07 
1,07 


407 
407 
1,03 


403 
£03 
L102 


101 
1,00 
400 


399 
398 
398 


397 
396 
395 


a> 
be 
391 


390 
a0 
389 


(2) 
OP 
OR (FICE 
“> O 


091 
099 
1.30 


2409 
3697 
10.1 


12.7 
13.9 
Lb 


17 1 
18.8 
20.3 


24.0 
2709 
3242 


36.5 
41.2 
450d 


51.4 
2309 
D9 oD 


63.7 
68.5 
76.0 


520 
520 
§20 
520 
520 
520 
§20 
520 
520 
520 
520 
520 
52h, 
525 
525 
525 
530 
233 


235 


(4) 
Mea 


ACTUAL 
£8 


SEC 


066 
ofl 
82 


1.03 
1.41 
Lies 


ZelyD 
2.60 
2e/l 


(5) 
V 


ET. 
EC 


20.0 
20.7 
ao eo 


30,2 
41.4 
G5 «1 


7320 
763 
7965 


$39 
87.6 
91.0 


995 


106.8 
113.3 


121.0 
12722 
133.0 


140.7 
143.9 
151.0 


156.0 
161.2 
16542 


(6) 
0] 


e i, 7 
o51 
268 


(7) 
AP 


METER 


“HO 


FRIMARY 


032 
o32 
039 





Barometer 29.85 "he 


(8) 
V/ 


RPM 











Ce 
atl 


¥ 
3 


AN] 
eo 8 @ 


® 
Wwnt WwOnWg GONM LTE" O WUOsF 


Wr Un A Vr A An 
Oe ™a CN 
i] t 


Lr nu Wn WIN 
e e 


SINE NOD CONN 


TAA AR 


(30) 


ie. 
DO xy 
iil 


PRIMARY 


07 
097 
lagi 


Lids. 
‘ep 

051 
-~ ela 
= 62 


-./°? 


-1.20 
aoe 5°) 
=i é 56 


=.96 
037 
1.80) 


2 otitis 
Br oe 
bye ly? 
be 52 
bo 35 
i oe 


led 5 
tie > 
3.76 


301 
3eke 


TABLE I 


(continued ) 
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2 AQ ww) SVvid 
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~el5 
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aa} 
0 iy 


(13) 
K 


Ae ol 
a 
43h 


4.32 
5.37 
28.1 


~13i, 
~27.0 


m2lne3 


a15 e ly 
=e e x 
-12 e § 


“2265 
62.8 
13a 


10.8 
7o3n 


Or hO 







(14) (75 
AP daa 
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(1) 
FP ABS 
ORIFICE 
“HO 


(2) 
aP 
ORIFICE 
“te oO 


099 
1.10 
1.46 


2036 
3023 
Le 
780 
12.25 
1/,.30 


3627 
me: 
50. 9 


567 
67.5 


FP 2 
90.6 


(4) 
Ma, 


ACTVAL 
<8 


SEC 


0/0 
o7 
085 


1.05 
1.26 
Le52 


1.95 
co4l 
2660 


20/0 
285 
2096 


3.10 
320 
3033 


3259 
3.83 
4.06 


4.31 
1 oh5 
be 69 


eI 
5428 
5.67 
5.99 


TABLE II 
Flovmeter Data For 4350 RPM 


(5) 
V 


Bec 


21.05 


2202 
20 ke 


3243 
37 07 
L507 


53 «3 
7202 
779 
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(7) 

LBP 
METER 

“tho 
PRIMARY 


47 
oh f 
055 


072 
083 
e931 


086 
064 
032 


045 
e3l 
030 





Barometer 29.76 "hg 


(8) 
N 


R PWM 
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LP, n# 
Ww x/O 
PRIMARY 


1.08 
1.11 
1.27 


1.66 
1.90 
2009 


1.98 
1.47 
1.20 


1.04 
0/2 
069 


046 
028 
oli 


028 
1.46 
2 TS 


lye 2k, 
Dend 


6635 


702 
7.80 
3659 


8.97 


(11) 
xz 
V 


3.60 
+ Le 
2097 


203) 
2.02 
1.66 


1.30 
16> 
097 


93 
685 
~82 
79 
76 


069 
¢66 
062 


060 
Sek 
oDk 


050 
048 
oy. 


o2 
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(12) 
ae 


& 


4063 
4.28 
3.79 


302 
2055 
1.90 


1.11 
oD 
236 


030 
019 
o17 


10 
206 
03 


005 
021 
03D 


051 
056 
61 


060 
059 
056 


03 


It 


(13) 
K 


(14) ies 
AP 4P x) 
ete Wee 
“H,O0 SECONDARY 
SECONDARY 

° 1.72 
075 1.72 
o87 2-00 
1.10 2053 
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APPENDIX II 





ORIFICE DATA 


Actual mass rate of flo’: for use as a standard by 


which to judge the flowmeter is computed fron the follow- 





ing equation: * 


Ma = HFS D, KY /S DP RIFICE (7 


where M,= mass rate of flow, 1lbm/sec 


Do= orifice diameter = 6.30 inches 
K = flow coefficient, cimensionless 
¥Y = expansion factor, dimensionless 


P,= static pressure before orifice, in. hg. abs. 
T = temperature before orifice, °R 


MAP orifice = pressure drop across orifice, in. water 


Py» T, andOP orifice are measured quantities 
recorded for each run at each mass flo increnent. 
Y, the expansion factor, is deternined fro Fig. 


IIeA which is a rlot cof the relation # 


=e =1-/.4/ va 35(B) }f Bese] (II~+2) 


me leery, Ref. 15, p. 9, Eqs. 





# Leary, Ref. =, De aye ergs 
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which, with constant diameter ratio, eS, and constant 
( 


ratio of specific heats, k, becomes a linear function of 
ADPoriric € QD Poriric € 


» Therefore, with a value of con~- 
5 A 
puted from measured quantities, Y may be read from Fig. 
ITA. 


K, the flow coefficient, has been experimentally 
determined by the ASME for this type of orifice and has 
been plotted versus Reynolds Number for specific values 
of 5s * . Use of a plot such as this would require a 
trial and error computation of K since the mass flow 
cannot be found until K is known and X cannot be found 
until the Reynolds Number, which is a function of mass 
flow, is known. 

To eliminate the necessity for a trial and error 
solution for K, two assumptions have been made: 

(a) Py, assumed constant = 29.9 in. hg. = 407 in 

water 

(b) T assumed constant = 530 °R 

With these assumptions Y becomes a function of 


4OPorifice only and equatioyg (II-2) becomes 
Y= /-, 000798 APoririce (II~3} 





* Ref. 25, Fig. 34¢ 
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Substituting this relation for Y into equation (II-1) 


we obtain 


Ma, = (4) 6.3) CKN7-. 200798 APogitce) |S [APoritice - 
which reduces to 


II- 
M, =/ 08k (/-. 000798 APoesrice)/APoririce (II-4) 





The Reynolds Number based on the orifice diameter 


is given by * 
a _ 15.28 Ma (II-5) 
< De 
Substituting equation (II-4) and taking “x , the vis- 
cosity = .000012 lb/ftesec, the expression for Reynolds 


Number becomes 


he =(2 /6 Ko) K (/-. 0007 98 AReairice) (AP ow i1ce Dae 


For a direct method of finding K, a curve of K 
versus OPorifice Would be desirable. Therefore, we pro- 
ceed as follows: 

(a) For an assumed value of OP. psrice, assume a 

K and compute R, from equation II-6. 
(b) Read K for computed Re from Fig. 34f of Ref. 





* Leary, Ref. 15, pe 12, Eqe 10d 









25 and use as new assumed K for repeat of 





step (a). 
(c) Plot K versus APoyiricg for various assumed 


values of 4 Porificee This plot will be a 









direct means for determining K, the flow 


coefficient. 





Firpure II-B is the result of the above procedure. 
Test calculations using Fig. 34f of Ref. 25 exclusively 


show that, for the range of pressures and temreratures 







investigated, there is negligible error introduced by 
assuming a constant P, anc T for the purpose of deter- 
mining «. 

With Y and K thus detemnined from Figs. II-A and 
II-B, actual mass rate of flow is computed from equation 


(II-1). 
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Fig. 2a Uniform Potential Flow Around 
A Circular Cylinder 


Fig. 2b Potential Vortex Flow 
With Cylinder As Core 


Fig. 2c Superposition Of Two 
Preceding Flows 
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From L. Prandtl] and O. Tietjens, Hydro. and Acromechanik, 2 (Springer, Berlin) 


VYV-CYL FANG SPEED 
UYU: STREAM VEX. 


b 


From L. Prandtl and 0. Tietjens, Hydro. und Acromechanik, 2 (Springer, Berlin) 
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